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ABSTRACT 

Nearby dwarf galaxies exhibit tight correlations between their global stellar and dynamical prop- 
erties, such as circular velocity, mass-to-light ratio, stellar mass, surface brightness, and metallicity. 
Such correlations have often been attributed to gas or metal-rich outflows driven by supernova energy 
feedback to the interstellar medium. We use high-resolution cosmological simulations of high-redshift 
galaxies with and without energy feedback, as well as analytic modeling, to investigate whether the 
observed correlations can arise without supernova-driven outflows. We find that the simulated dwarf 
galaxies exhibit correlations similar to those observed as early as z ~ 10, regardless of whether super- 
nova feedback is included. We also show that the correlations can be well reproduced by our analytic 
model that accounts for realistic gas inflow but assumes no outflows, and star formation rate obeying 
the Kennicutt-Schmidt law with a critical density threshold. We argue that correlations in simulated 
galaxies arise due to the increasingly inefflcient conversion of gas into stars in low-mass dwarf galax- 
ies rather than supernova-driven outflows. We also show that the decrease of the observed effective 
yield in low-mass objects, often used as an indicator of gas and metal outflows, can be reasonably 
reproduced in our simulations without outflows. We show that this trend can arise if a significant 
fraction of metals in small galaxies is spread to the outer regions of the halo outside the stellar extent 
via mixing. In this case the effective yield can be significantly underestimated if only metals within 
the stellar radius are taken into account. Measurements of gas metallicity in the outskirts of gaseous 
disks of dwarfs would thus provide a key test of such explanation. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - galaxies: abundances 
- galaxies: fundamental parameters ~ galaxies: dwarf 



1. INTRODUCTION 

Observational studies of galaxies have revealed the ex- 
istence of tight correlations between their global stellar 
and dynamical properties, such as circular velocity, mass- 
to-light ratio, stellar mass, surface brightness, and metal- 
licity. For example, scalings of the mean metallicity of 
a galaxy with other global properties such as its lumi- 
nosity, stellar mass, or total mass, have been established 
for galaxies of a large range of masses and morphologies 
(e.g., Lequeux et al. 1979; Garnett & Shields 1987; Zarit- 
sky et al. 1994; Garnett 2002; Prada & Burkert 2002; 
Dekel & Woo 2003; Tremonti et al. 2004; Pilyugin et al. 
2004; Simon et al. 2006; Lee et al. 2006a). Theoretical 
modeling of these correlations can help us to identify key 
physical processes shaping global properties of galaxies. 

Supernova energy feedback to the interstellar medium 
(ISM) with putative associated gas outfiows from galax- 
ies has long been a favored mechanism to explain these 
correlations and other properties of low-mass galaxies 
(e.g., Larson 1974; Dekel & Silk 1986; Arimoto & Yoshu 
1987). The strong effect of the SN feedback on proper- 
ties of low-mass galaxies is also a standard assumption 
of the semi-analytic models of galaxy formation (e.g., 
Lacey et al. 1993; Kauffmann et al. 1993; Cole et al. 
1994; Somerville & Primack 1999; Benson et al. 2003; 
Croton et al. 2006). Dekel & Woo (2003) argue, for 
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example, that the correlations between circular veloc- 
ity, stellar mass, metallicity, and stellar surface densities 
they find for the nearby dwarf galaxies can be repro- 
duced with a simple semi-analytic model incorporating 
the effect of supernova feedback on the gas component 
of dwarf galaxies. 

The existence of galactic winds has indeed been obser- 
vationally established in star-forming galaxies at low and 
high redshifts (e.g., Heckman et al. 2000; Pettini et al. 
2001; Martin et al. 2002; Strickland et al. 2004; Martin 
2005; Ott et al. 2005, see Veflleux et al. 2005 for a recent 
review). An indirect evidence for metal loss in winds is 
the high-metallicity of the diffuse intergalactic medium in 
groups and clusters (e.g., Renzini et al. 1993). However, 
the extent to which the outflows affect the global prop- 
erties of galaxies and whether the outflow gas escapes 
the gravitational potential of its host halo or rains back 
down to the disk remains uncertain. Mac Low & Ferrara 
(1999) and D'Ercole & Brighenti (1999) used numerical 
simulations to show that global gas blowout is inefficient 
in all but the smallest mass halos [M < 10'' Mq; see also 
Ferrara & Tolstoy 2000; Marcolini et al. 2006), although 
metal-enriched SN ejecta may be removed efficiently. It 
is also uncertain whether correlations of metallicity with 
galaxy stellar mass and mass-to-light ratios can be at- 
tributed solely to such winds. 

In the context of cosmological simulations, the SN feed- 
back is usually ineffective, unless ad hoc phenomeno- 
logical recipes enhancing SN feedback are employed 
(e.g., Navarro & White 1993; Navarro & Steinmetz 
1997; Thacker & Couchman 2001; Marri & White 2003; 
Springel & Hernquist 2003; Scannapieco et al. 2006). It 
is likely that this inefficiency is at least partially due 
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to the inability of current cosmological simulations to 
resolve the scales and processes relevant to stellar feed- 
back. However, it is also possible that the actual effects 
of feedback on the global properties of galaxies are fairly 
small in reality. Cosmological simulations at this point 
cannot make ab initio predictions about the importance 
or inefficiency of stellar feedback. Observations, on the 
other hand, although showing evidence for the presence 
of large-scale winds in actively starforming galaxies (e.g., 
see Veilleux et al. 2005, for review), are rather uncertain 
in their estimates of wind mass loss to provide reliable 
direct obscrvationally-motivatcd feedback recipes for im- 
plementation in simulations. Nevertheless, we can gauge 
the importance of supernova energy injection by compar- 
ing galaxies formed in simulations performed with differ- 
ent assumptions about feedback to observations. In this 
respect, the smallest dwarf galaxies present the ideal test 
case because they can be expected to be the most sus- 
ceptible to the effects of SN feedback due to their shallow 
potential wells (Larson 1974; Dekel & Silk 1986). 

The role of feedback in the formation of dwarf galaxies 
has been investigated in several recent stiidics. Tassis 
et al. (2003) used Eulerian AMR simulations of galaxy 
formation with star formation and SN feedback and 
found that correlations between mass-to-light ratio and 
metallicity and stellar mass and metallicity, similar to 
those observed for the nearby dwarfs, arc exhibited by 
dwarf galaxies in their simulations at z >3. The simula- 
tions of Tassis et al. (2003) used a rather extreme amount 
of energy per supernovae to maximize the effects of feed- 
back. They did not test, however, whether feedback or 
some other mechanism is in fact the dominant factor 
in shaping the correlations. More recently, (Kobayashi 
et al. 2006) found a tight relation between stellar metal- 
licity and stellar mass at all redshifts in their Smooth 
Particle Hydrodynamics (SPH) cosmological galaxy for- 
mation simulations, which they attributed to the mass- 
dependent galactic winds. The winds in the model of 
Kobayashi et al. (2006) do affect the gas and metal con- 
tent of galaxies significantly, with the effect increasing 
towards lower mass systems (see their Fig. 16), which 
plays an important role in shaping the resulting correla- 
tion of metallicity and stellar mass. 

De Rossi et al. (2006) have recently studied the ori- 
gin of luminosity- and stellar mass-metallicity relations 
in cosmological simulations. They find that correlations 
similar to those observed arise already at high redshifts, 
although supernova feedback in their simulations is not 
efficient. This is consistent with our results presented be- 
low. Simon et al. (2006) seeked observational correlations 
between metallicity, surface brightness and mass-to-light 
ratio within a single galaxy (M33) and found that in- 
deed such correlations exist. They used SPH simulations 
and analytical calculations to show that feedback is not 
responsible for these correlations. They found that the 
presence of feedback decreases the scatter in such rela- 
tions, although this may at least in part be due to their 
numerical implementation of feedback. Most recently. 
Brooks et al. (2006) presented a study of the origin and 
evolution of the mass-metallicity relation in cosmological 
simulations. They also concluded that the relation arises 
not due to the mass loss in winds, but due to increasing 
ineflaciency of star formation in smaller mass galaxies. 
However, they argue that low effective yields of dwarf 



galaxies can only be explained by the loss of metals in 
winds. We argue below that there is an alternative ex- 
planation for the low effective yields of dwarfs. 

In this paper we investigate the origin of the observed 
mass-metallicity and other correlations of dwarf galax- 
ies with the specific goal of testing the role of SN feed- 
back. In particular, we discuss results of several high- 
resolution cosmological simulations of galaxy formation 
started from the same initial conditions and run with 
the same prescriptions for star formation and metal en- 
richment of the ISM, but with different models for gas 
cooling, UV heating (optically thin approximation vs. 
self-consistent radiative transfer of the UV radiation), 
and supernova feedback (see § 4 for details). We follow 
the simulations until z ^ 3. 

Although our models simulate galaxies at high red- 
shifts, we study the evolution of the correlations with 
redshift from z = 9 to z ~ 3 and show that the trends 
with metallicity and stellar mass arc approximately pre- 
served as the galaxy population evolves allowing extrapo- 
lation of the results to the present epoch. Therefore, our 
resiilts and conclusions give us insight to the processes 
responsible for the properties of observed low-redshift 
dwarfs. There is indeed observational indication that 
the stellar mass-metallicity correlation, for example, is 
already established by 2; ~ 1 — 2 (Kobulnicky & Kew- 
ley 2004; Savagho et al. 2005; Erb et al. 2006), although 
the metallicities of galaxies of a given stellar mass ap- 
pear to be lower by w 0.3 dex at z ^ 2 compared to 
z = (Savaglio et al. 2005; Erb et al. 2006). In addition, 
for some of the Local Group dwarf galaxies that exhibit 
these correlations, a significant fraction of their stellar 
mass was built at z > 3 (e.g.. Dolphin et al. 2005). The 
correlations for these dwarfs should thus be largely set 
at high redshifts, which we probe with our simulations. 

Our main result is that correlations similar to those ob- 
served can be reproduced in all runs, which suggests that 
supernova energy feedback is not required to explain the 
observed properties of dwarfs. To interpret the simula- 
tion results we use a simple open-box analytic model for 
the evolution of the baryonic component of dwarf galax- 
ies. We show that a model with gas accretion but no 
outflows and the Kennicutt law for star formation with 
the critical density threshold for star formation repro- 
duces the results of our simulations. 

This paper is organized as follows. In §2 we summa- 
rize the observationally established correlations between 
global properties of dwarf galaxies that we wish to inter- 
pret. Our analytic model and chemical evolution calcu- 
lations are described in §3. In §4 we describe in detail 
our cosmological simulations, and in §5 we discuss how 
global properties of the simulated galaxies are correlated 
with each other, and we compare these correlations to 
the observed ones, and to the analytic model of §3. A 
discussion of observational results on the effective yield 
and corresponding constraints on the existence of out- 
flows is given in §5.4. We discuss our findings in §6, and 
we summarize our conclusions in §7. 

2. OBSERVED CORRELATIONS 

Correlations of observables in dwarf galaxies relate the 

stellar mass, M,, with the metallicity Z, the surface 
brightness /x* and the maximum circular velocity Vm of 
dwarf galaxies. These can be parametrized as the power 
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law relations: 

Z(xM:^; /X, ocMr"; (x M^^ . (1) 

The existence of these correlations (in several cases in 
the form of correlations between absolute magnitudes, 
rather than stc;llar masses, and Z. fi^ and ^1,1) li^is been 
established through observations of galaxies with a large 
range of masses and in a variety of settings, including 
nearby dwarfs (e.g. Skillman et al. 1989; Dekel & Woo 
2003 based on observations taken from Mateo 1998, van 
den Bosch 2000 and references therein; Lee et al. 2006a; 
van Zee & Haynes 2006), dwarf galaxies in the Sloan Dig- 
ital Sky Survey (SDSS; Tremonti et al. 2004; Kauffmann 
et al. 2003; Blanton et al. 2003; Bernard! et al. 2003), the 
Hubble Deep Field (Driver 1999), the Ursa Major cluster 
(Bell & de Jong 2001) and the Virgo cluster (Ferguson & 
Binggeli 1994). However, the exact values of the slopes 
of these relations are affected by systematic uncertainties 
involved in obtaining stellar masses from the observed 
quantities (magnitudes). In general, the slopes depend 
on the the mass range in which they are determined. 

For example, in the case of the Z — scaling, Dekel 
& Woo (2003) find that for the Local Group dwarfs in 
the mass range of lO^M© < M* < 1O^°M0, nz ~ 0.4, 
while the SDSS data of Tremonti et al. (2004) (see also 
Gallazzi et al. 2005) indicate that nz increases with de- 
creasing M,: it is quite flat for > 3 x lO^^M©, while 
it appreciably exceeds the Dekel & Woo (2003) value by 
M ~ few X lO^M©. Most recently, Lee et al. (2006a) find 
nz ~ 0.29 ± 0.03 for the dwarf galaxies in their sample. 

In the case of the /i* — M* scaling, Dekel & Woo (2003) 
find for the Local Group dwarfs the slope of ~ 0.55, 
although their data points exhibit a scatter around the 
best-fit line which is large compared to their adopted er- 
ror bars. For the same scaling, Kauffmann et al. (2003) 
using SDSS data for the surface density of stellar mass 
within the half-light radius, find that for a given value of 
there is appreciable scatter (about an order of mag- 
nitude) in ;Li* , and the median value of /i* scales roughly 
as for 10« < < 3 x lO^^Mo, and it becomes 

almost independent of M* for larger M,. 

In the case of the M* — Vm scaling (the TuUy-Fisher 
relation). Bell & de Jong (2001) using galaxies with 
3 X IO^Mq < < 1O"M0 find ny ~ 0.21 - 0.24, and 
a shallower correlation with baryonic mass (i.e. stars 
plus gas), ny ^ 0.27 — 0.30. Several recent studies 
also find that the Vm — M^, relation for low-mass galax- 
ies (T4i £ 90 km/s) is steeper and exhibits considerably 
more scatter than the baryonic Tully-Fisher relation of 
the same galaxies (Matthews et al. 1998; McGaugh et al. 
2000; Bell & de Jong 2001; Gurovich et al. 2004; Mc- 
Gaugh 2005; Geha et al. 2006). These studies find slopes 
of « 0.1 — 0.2 for the M» — relation of the dwarf 
galaxies (McGaugh et al. 2000; Gurovich et al. 2004) and 
nv ~ 0.25 — 0.33 for the baryonic Mbaryon — rela- 
tion (e.g., McGaugh et al. 2000; Verheijen 2001; Bell & 
de Jong 2001; McGaugh 2005; Geha et al. 2006). For 
the Local Group dwarfs, Dekel & Woo (2003) obtain 
nv ^ 0.37 for lO^M© < M* < IO^^Uq and ny ~ 
for IO^Mq < < IO^Mq. 

Finally, another tight relation found in Local Group 
dwarf galaxies is the correlation between the dynami- 
cal mass-to-light ratio and average metallicity (Prada & 



Burkert 2002): 

log(Mdy„/Ly)oc-[Fe/H]. (2) 

3. ANALYTICAL MODEL FOR THE EVOLUTION OF THE 
BARYONIC COMPONENT OF DWARF GALAXIES 

In this section we present a simple analytic model for 
the evolution of the baryonic component (gas, stars and 
metals) of dwarf galaxies at high redshifts, which we will 
make to interpret our simulation results. The model is 
based on the following assumptions: 

(1) Galaxies accrete gas at a rate consistent with sim- 
ulations, but no gas or metals escape outside the virial 
radius of the parent halo. 

(2) The gas mass of each galaxy is proportional to its 
total (dynamical) mass, Mgas oc M . At early epochs 
M^, <C Mgas, which implies for the baryonic mass Mb = 
Mgas + M* K, Mgas- The gas fraction within halos in 
our simulations is approximately independent of the to- 
tal (dynamical) mass for M^oi ^ lO^'M© (see §5, Fig. 
2), which is consistent with the results of independent 
simulations of Tassis et al. (2003). Therefore this as- 
sumption is an excellent approximation for all but the 
smallest objects in the simulation, which are excluded in 
our analysis of other global quantities. Although at lower 
redshifts the total baryon fraction in small mass galaxies 
M < 10^" Mq can be suppressed due to heating by the 
cosmic UV background (e.g., Gnedin 2000; Hoeft et al. 
2006), our results would still be applicable to galaxies in 
which a significant fraction of stars is formed at high red- 
shifts before such suppression has occurred in the mass 
range we consider. We should stress that we intend to 
apply our model only to such stellar populations. 

(3) The gas in galaxies forms an exponential disk with 
radial extent proportional to the virial radius of the par- 
ent halo. The disk is assumed to be more extended and 
consequently less compressed for smaller objects with 
Tvir < 10** K, in which gas cannot cool efficiently. 

(4) Gas is converted into stars at a rate given by the 
Kennicutt law of star formation with a density threshold 
(e.g., Martin & Kennicutt 2001). Gas at densities below 
the threshold does not form stars. Note that the Ken- 
nicutt law with the threshold is reproduced by our nu- 
merical simulations, but it is not assumed a priori. The 
numerical simulations implement a much simpler recipe 
for star formation (see §4). 

(5) Instantaneous recycling approximation: i.e., no 
time delay is assumed between the birth of a genera- 
tion of stars and the corresponding metal enrichment of 
the ISM. Hence, we implicitly assume that the metal en- 
richment due to type II supernovae is dominant. This is 
appropriate for all but the most massive dwarf galaxies. 

The key ingredient of our model is the relation between 
the gas mass and the stellar mass of a galaxy. We will 
first derive this correlation, and subsequently use it to 
derive relations between the stellar mass, gas mass, stel- 
lar surface density, metallicity, and circular velocity of 
high-z galaxies. 

3.1. Stellar Mass - Gas Mass Relation 

The relation between stellar mass and total gas mass 

in dwarf galaxies is central to understanding the origin of 
correlations between other observable quantities'"^. In this 

^ This relation may not hold for some galaxies, if, for instance, 
their gas is lost due to tidal or ram pressure stripping. 
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Fic. 1.- Sensitivity of the analytic model to the angular mo- 
mentum parameter A and the surface density threshold for star 
formation, S^i^, at z = 4. Solid line: fiducial model with A = 0.033 
(its most probable value) and Stjj = 5Mqpc~^. Dashed lines: 
model results for Sth = 5Mqpc^^ and A = 0.023 (left dashed 
line) or A = 0.044 (right dashed line). Dotted lines: model re- 
sults for A = 0.033 and Sth = 2M0pc-2 (left dotted line) or 
Sth = 10 Mq pc~^ (right dotted line). 



section, wc derive stich a relation. In subsequent sections, 
we will use it to explain the correlations discussed in §2 
and we will directly compare it against our simulation 
results. 

The rate of change of the gas mass, assuming that there 
are no significant gas outflows due to SN-driven winds, 
is 

^(Mga,) = F-aV, (3) 

where F is the gas mass accretion rate, ijj is the star 
formation rate, and a is the fraction of stellar mass locked 
in long-lived stars and remnants (i.e., not returned to the 
ISM). The rate of change of the stellar mass is 



dt 



(4) 



.^'-^gas, with a slope set by the slope of the star- formation 
law. At lower Mgas, this relation exhibits a break, which 
is a result of the threshold in the star-formation law re- 
sulting in very inefficient star formation at the low den- 
sities achieved by inefficiently cooling low-mass galaxies. 

Remarkably, we find that the best fit to our data at 
= 4 is obtained when both free parameters have their 
fiducial values, A = 0.033 (the most probable value) and 
Sth = SMqPC"^ (the "canonical" observational value). 
Figure 1 shows the sensitivity of our analytic model to 
the values of A and Sth- The solid line represents the z = 
4 model results for the fiducial parameter values (which 
are also the best- fit parameters). The dashed lines show 
the range of A values that brackets the scatter of our 
simulation data points at 2 = 4. The range of A values 
represented here is 0.023 ^ A < 0.044 (which contains 
~ 40% of the area tmder the A— probability distribution, 
and is centered aromid the distribution peak). In this 
case, Sth is kept at its fiducial value. The dotted lines 
show the range of Sth values that brackets the scatter 
of our data. The range of Sth values represented here is 



2M 



0PC 



-2 < 



Sth ^ IOMqpc^ . In this case, A is kept 
at its fiducial value. The asymptotic slope of the curve at 
higher masses, remains largely independent of parameter 
choice. 

3.2. Mass-to-light ratio vs. metallicity 

Under the assumptions of negligible infiow or outflow 
of metals and instantaneous recycling, and for a fixed 
initial mass function, the mass of metals of \ow-Z, high- 
z galaxies is proportional to the stellar mass^. From 
Mz = ZMgg,s oc M* we can immediately write 



Z oc 



(5) 



gas 



If the dependence of tp and F on Mgas, as well as any 

explicit dependence on redshift, are known, then equa- 
tions (3) and (4) can be integrated in time as a system 
of two ordinary differential equations. Different initial 

conditions for Mgas at some early time (for which we can in combination with Eq. (5) and assumptions Mgas oc M 



Assumption (2), Mgas M, then gives Z cx M^/M. 
The luminosity of an object, L, can be considered pro- 
portional to its stellar mass^. Then, Eq. (5) becomes 
Z (X L/M or 

M/LocZ-^, (6) 

in agreement with observations and, as shown in §5, the 
results of our simulations. This is a simple interpretation 
of the correlation between the mass-to-light ratio and 
metallicity found by Prada & Burkert (2002). 

3.3. Additional Correlations 
We can use the M„as — M* relation we derived in 53.1 



safely assume that M* = 0) will produce different pairs 

of (Mgas,M») at any time of interest. In this way, we 
calculate parametrically the prediction of our analytic 
model for the M* — Mgas relation. 

As we show in the Appendix, F is proportional to the 
gas mass and V depends on the extent of the gaseous 
disk (which in turn depends on the dynamical mass of 
the galaxy and the angular momentum parameter, A) 
and the star formation density threshold, Sth- 

The free parameters which control the normalization 
and the shape of the Mgas — relation are A and Sth- 
The solution of Eqs. (3) and (4) for the fiducial values of 
these parameters is plotted with the solid line in Fig. 1. 
At high gas masses, M* has a power-law dependence on 



and L oc M^ to derive additional relations between stellar 
mass and stellar surface density, metallicity, and circular 
velocity of high-z galaxies. We will overplot these re- 
lations with our simulation results and discuss them in 
§5- 

^ For this proportionality to hold, Z needs to be low enough so 
that the mass of freshly produced metals is much laxger than the 
mass of metals locked in long-lived stars and stellar remnants. This 
can be shown to be true for Z <C Zq for the metallicity enrichment 
implemented in our simulations, and which always holds in our 
simulated galaxies. 

This approximation is appropriate for older stellar populations. 
Wc use it here because we assume that the stellar populations in 
our high-z galaxies would actually be observed today, when they 
have aged substantially. 
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TABLE 1 
Simulations 



Simulation* SN energy feedback cooling 3-D rad. transfer 



FEC yes equilibrium no 

NFEC no equilibrium no 

FNEC-RT yes non-equilibrium yes 

F2NEC-RT yes non-equilibrium yes 



Throughout the text we denote different runs with abbreviations indicat- 
ing the different implementations of key physical processes. In the abbrevi- 
ations: F is for feedback, NF is for simulations with SN feedback explicitly 
turned off, and F2 is for an alternative implementation of feedback with de- 
lay after the star formation event; EC is for equilibrium cooling using Cloudy 
tables, and NEC is for non-equilibrium cooling and reaction network of ionic 
species of H, He, and H2; RT indicates runs with self-consistent radiative 
transfer. 



4. NUMERICAL SIMULATIONS 

In this study we use four simulations (Table 1) of the 

early {z > 3) stages of evolution for a Lagrangian re- 
gion of a MW-sized system of total (dark matter and 
baryons) virial mass « 10^^ Mq at z = 0. The 
simulations were performed using the Eulerian, gasdy- 
namics -t-AT-body Adaptive Refinement Tree (ART) code 
(Kravtsov et al. 1997; Kravtsov 1999). In both the gas- 
dynamics and gravity calculations, a large dynamic range 
is achieved through the use of adaptive incsh refine- 
ment (AMR). At the analyzed epochs, the galaxy has 
already built up a significant portion of its final total 
mass: 1.3 x lO^^/i"^ at z = 9 and 2 x lO^'^h'^ Mq at 
z = 4. The evolution is started from a random realiza- 
tion of a gaussian density field at z = 50 in a periodic box 
of 6h~^ Mpc with an appropriate power spectrum and is 
followed assuming flat ACDM model: Qq = 1 — f^A = 0.3, 
fib = 0.043, h = iJo/lOO = 0.7, n, = 1, and erg = 0.9. 

To achieve the high mass resolution, a low-resolution 
simulation was run first and a galactic-mass halo was se- 
lected. A Lagrangian region corresponding to five virial 
radii of the object at z = 0, corresponding to a re- 
gion of ^ "ihr^ comoving Mpc in diameter, was then 
identified at z = 50 and re-sampled with additional 
small-scale waves (Navarro & White 1994, the specific 
procedure described here is described in Klypin et al. 
2001). The total number of DM particles in the high- 
resolution Lagrangian region is 2.64 x 10® and their mass 
is mDM = 9.18 X 10^/i~-^ Mq. In addition to the main 
progenitor of the MW-sized system this Lagrangian re- 
gion contains several dozens of smaller galaxies spanning 
a wide range of masses. We will use these galaxies to 
analyze correlations between their properties. 

The code used a uniform 64^ grid to cover the en- 
tire computational box. The Lagrangian region, how- 
ever, was always unconditionally refined to the third re- 
finement level, corresponding to the effective grid size of 
512^. As the matter distribution evolves, the code adap- 
tively and recursively refines the mesh in high-density 
regions beyond the third level using two refinement cri- 
teria: gas and DM mass in each cell. A mesh cell 
was tagged for refinement if its gas or DM mass ex- 
ceeded 1.2 X W'hr^ Mq and 3.7 x W'hr^ Mq, re- 
spectively. The maximum allowed refinement level was 
'max = 9. The volume of high-density cold star form- 
ing disks forming in DM halos was refined to 'max = 9. 



The physical size of mesh cells in the simulations was 
Axi = 26.16 [10/(1-Fz)]29-' pc, where I is the cell's level 
of refinement. Each refinement level was integrated with 
its own time step Af; = Ato2-' w 2 x 10^ 2^"' yrs. The 
time steps are set by a global Courant-Priedrich-Levy 
condition. 

The analyzed simulations were started from the same 
initial conditions, but evolved using different assump- 
tions about cooling and heating processes accompanying 
galaxy formation. All four runs include star formation 
using the recipe described in Kravtsov (2003). Namely, 
the gas is converted into stars on a characteristic gas 
consumption time scale, r*: p* = Pg/r*. We use con- 
stant, density- independent r*. This is different from the 
commonly used density-dependent efficiency, but may be 
more appropriate for star forming regions on ^ 100 pc 
scales which arc; resolved in the simulations where con- 
stant star formation efficiency is indicated by observa- 
tions (e.g.. Young et al. 1996; Wong & Blitz 2002). As 
shown by Kravtsov (2003), such constant efficiency as- 
sumption at the scales of molecular clouds results in the 
Kennicutt-likc star formation law on kiloparsec scales, 
although a Kennicutt law is not a priori assumed. The 
star formation was allowed to take place only in the 
densest cold regions, pg > psF and Tg < Tsf, but no 
other criteria (like the collapse condition V • v < 0) 
were imposed. We used t» = 4 Gyrs, Tsf ~ 10** K, and 
PSF = 1-64 Mqpc~^ corresponding to atomic hydrogen 
number density of nH,SF = 50 cm^"^. 

Our fiducial simulation (FEC) includes metallicity- 
dependent cooling and UV heating due to cosmological 
ionizing background with equilibrium cooling and heat- 
ing rates tabulated for the temperature range 10^ < T < 
10'^ K using the Cloudy code (ver. 96b4, Ferland et al. 
1998). The cooling and heating rates take into account 
Compton heating/cooling of plasma, UV heating, atomic 
and molecular cooling. It also includes energy feedback 
and chemical enrichment due to supernovae assuming the 
Miller & Scalo (1979) stellar initial mass function (IMF) 
and stellar masses in the range 0.1 — 100 Mq. All stars 
with M* > 8 Mq deposit 2 x 10^^ ergs of thermal energy 
and a mass fzM^ of heavy elements in their parent cell 
(no delay of cooling was introduced in these cells). The 
metal fraction is /z = min(0.2, O.OlAf* — 0.06), which 
crudely approximates the results of Woosley & Weaver 
(1995). 

The second simulation, NFEC, is identical to the FEC 
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run in all respects, except that it did not include the en- 
ergy injection due to SNe (chemical enrichment due to 
SNe is still included). Note that the cooling rates in the 
run with feedback accounted for the local metallicity of 
the gas, while in the run with no feedback the signifi- 
cantly lower zero- metallicity cooling rates were used®. 

The third simulation. FNEC-RT, includes star forma- 
tion and SN enrichment and energy feedback in the same 
way as the FEC run, but uses self-consistent 3-D ra- 
diative transfer of UV radiation from individual stellar 
particles using the OTVET algorithm (Gnedin & Abel 
2001; Iliev et al. 2006) and follows non-equilibrium chem- 
ical network of hydrogen and helium species (the de- 
tails of the specific implementation of the OTVET algo- 
rithm on adaptively refined meshes will be described else- 
where). The simulation thus includes non-equilibrium 
cooling and heating rates which make use of the local 
abundance of atomic, molecular, and ionic species and 
UV intensity followed self-consistently during the course 
of the simulation. The metallicity-dependence of cool- 
ing rates is taken into account using optically thin equi- 
librium metal cooling functions from Sutherland & Do- 
pita (1993) in the high metallicity regime and from Pen- 
ston (1970) and Dalgarno & McCray (1972) in the low 
metallicity regime. The molecular hydrogen formation 
on dust is taken into account using Cazaux & Spaans 
(2004) rates, and gas cooling on dust is based on Draine 
(1981). 

A severe limitation of our current treatment of molec- 
ular hydrogen formation is that self-shielding of molecu- 
lar lines in the Lyman- Werner band is not taken into 
account, which results in a substantial underestimate 
of the molecular fractions at the high density regime 
(n > 100 cm~^). Note, however, that since our star for- 
mation criterion is based on the total rather than molec- 
ular gas density, this does not affect our star formation 
rates significantly. 

Finally, the fourth simulation, F2NEC-RT, is identical 
to the FNEC-RT but SN energy is injected with a delay of 
10'' years after the star forming event and energy release 
is spread over lO'' years after that. This variation is in- 
tended to test the sensitivity of our results to the specific 
details of the feedback implementation. Energy injection 
with delay was argued to enhance the efficiency and im- 
pact of the SN feedback on the ISM (Slyz et al. 2005), 
as the stellar particles have time to leave the densest re- 
gions where they form and release their energy at lower 
densities where the cooling time is longer. This change 
of feedback recipe results in a significant enhancement of 
the star formation efficiency. Perhaps surprisingly, the 
F2NEC-RT simulation forms nearly twice as much stellar 
mass as the simulation FNEC-RT. The reason is simple. 
With a delayed energy release the star forming regions 
are generally left intact and continuing forming stars in 
the simulation F2NEC-RT, while in the FNEC-RT the 
temperature of gas in star forming regions is increased 
due to SN feedback and star formation ceases for some 
period of time. 

We show below that the details of the specific imple- 

* This choice docs not appreciably affect the star-formation prop- 
erties of small galaxies, as shown by Kravtsov (2003), since, at least 
for the densities we are resolving here, the density probability dis- 
tribution function does not change significantly. 



mentation of feedback in our simulations do not affect the 
correlations of galaxy properties and our conclusions. 

5. RESULTS 

In this section, we discuss relations between global 
properties of galaxies in our cosmological simulations, 
how these depend on the detailed physics included in 
each simulation, and how they compare to the observa- 
tionally established correlations and our analytic models 
of §3. 

The panels in the figures discussed in this section are 
arranged in the following way: left panel - FEC simu- 
lation (SN feedback, uniform UV background, tabulated 
heating/cooling rates); middle panel - NFEC simulation 
(no SN energy feedback, uniform UV background, tab- 
ulated heating/cooling rates); right panel - FNEC-RT 
simulation (SN feedback, 3-D radiative transfer). We do 
not show the results of the F2NEC-RT simulations in all 
of the figures for clarity. In all cases, the results of the 
F2NEC-RT simulations are quite similar to those of the 
FNEC-RT run. 

5.1. Baryonic Content of Simulated Galaxies 

Figure 2 shows the gas and stellar mass fractions (to- 
tal gas mass over total dynamical mass, and stellar mass 
mass over total dynamical mass, respectively) within the 
virial radius as a function of the virial mass for all ob- 
jects in the simulations. For systems with M > 10^ M©, 
the gas fractions are approximately independent of mass, 
while at smaller masses a variety of effects (e.g., gas strip- 
ping, reionization) act to suppress the gas fractions and 
increase the scatter. This result is qiialitativcly con- 
sistent with the findings of the previous studies (e.g., 
Gnedin 2000; Chiu et al. 2001; Tassis et al. 2003; Hoeft 
et al. 2006). Figure 2 also shows that the stellar frac- 
tions for the redshifts and masses under consideration 
is an order of magnitude smaller than the gas fractions. 
This is due both to the inefficiency of star formation in 
these small systems and the fact that we analyze them 
at high redshifts, where they did not have the sufficient 
amount of time to convert their gas into stars. Most in- 
terestingly, the stellar fractions drop sharply in systems 
with M < 5 X 10^ M0. The objects with smaller masses 
do contain a universal share of the baryons in the form 
of gas, but fail to convert a sizeable fraction of it into 
stars. Note that all three simulations shown produced 
very similar results. This demonstrates explicitly that 
the inefficiency of star formation at small masses is not 
due to SN feedback. We note however that objects in 
these small-mass systems are resolved with just a few 
thousand DM particles. Convergence tests with much 
higher resolution simulations are therefore needed to test 
whether inefficient star formation at the smallest masses 
is not a numerical artefact. 

Figure 3 shows the star-to-baryonic mass ratio of the 
simulated galaxies plotted against their gas mass. A tight 
correlation between the two quantities exists, although it 
cannot be described by a single power law for the entire 
range of masses examined here. At the; high-mass cind, 
Mg > 10^ M0, the stellar mass is roughly proportional 
to the baryonic mass (and hence, in the gas-rich regime, 
to the gas mass). At lower masses, however, the rela- 
tion exhibits a break and becomes considerably steeper 
as the stellar mass of low-mass objects is suppressed (see 
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Fig. 2. — Gas and stellar mass fractions of objects in our simulations plotted against the total (dark matter + total baryonic) mass of 
each object. Left panel: simulation FEC; middle panel: simulation NFEC; right panel: simulation FNEC-RT. Open triangles show the 
gas fraction (Afgas/Aftot), while stars show the star fraction Af, /A/tot. Different colors correspond to different redshifts as detailed in the 
legend. The horizontal dashed line corresponds to the baryonic mass fraction {Q^^/Qq = 0.14) 




Fig. 3. — Star-to-baryonic mass ratio vs. gas mass for objects in each of the analyzed simulations. The gas mass here is the total gas 
in each halo. Left panel: simulation FEC; middle panel: simulation NFEC. Right panel: the FNEC-RT simulation. Different symbols 
correspond to different redshifts as detailed in the legend. The solid and dashed lines show the results of our analytic model for 2 = 4 and 
z = 9, respectively. 
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also Chiu ct al. 2001). This correlation is central in iden- 
tifying the origin of the scalings between the rest of the 
global properties we present below. The solid and dashed 
lines represent the prediction of the model discussed in 
§ 3, for redshifts z = 4 and z = 9, respectively. The 
agreement with the simulation data is excellent. Here 
again we see that the results of the three different simu- 
lations are consistent with each other. 

5.2. Correlations Involving Metallicities 

Before we proceed to discuss the correlation of the 
metallicities in the simulated galaxies with other observ- 
ables, we should make a note about the distinction be- 
tween stellar and gas metallicities. Metallicities in stars 
represent an average metallicity over the star formation 
history of the observed object, while gas metallicities are 
actually representative of the amount of metals currently 
present in the interstellar medium. Observationally, ei- 
ther one or both can be measured and used. For exam- 
ple, Dekel & Woo (2003) use stellar metallicities for dE 
galaxies and gas oxygen abundance for dirr, while Prada 
& Burkert (2002) use stellar metallicities in all cases. 
The metallicity in the gas and in the stars, although not 
the same in general, are related to each other. This is 
particularly true for the high-redshift galaxies in our sim- 
ulations. Figure 4 shows the metallicity in the gas versus 
the metallicity in stars for every object in our simulations 
at different redshifts. The solid line corresponds to the 
one-to-one relation Zgas = Z^,. Indeed, for the objects 
and redshifts we consider, the metallicity of stars tracks 
that of the gas very closely. 

A tight correlation observed for the nearby dwarf 
galaxies is the correlation of the dynamical mass-to-light 
ratio and metallicity (Prada & Burkert 2002). In § 3.2 we 
showed that such a correlation should naturally arise in 
low-metallicity dwarf galaxies without gas outflows and 
regardless of their star formation rate as long as their 
baryon mass is approximately proportional to their total 
mass. Figure 5 clearly shows that for galaxies in our sim- 
ulations the ratio of the total (baryonic + dark matter) 
to stellar mass is tightly correlated with metallicity, in- 
dependently of the details of star formation and cooling 
and inclusion of SN feedback. The form of the correla- 
tion, M/M* cx is similar to that exhibited by the 
observed dwarf galaxies. The correlation arises for ob- 
jects at the earliest epochs, which implies that a similar 
correlation should be expected to hold for the observed 
galaxies at high redshifts. 

Basic assumptions in the derivation of this correlation 
is that Mgas oc M and Z <^ Zq, implying that at the 
time this correlation is being established, the galaxy is 
in the gas-rich, low-metallicity regime. For galaxies in 
which star formation is inefficient and which have not 
lost substantial amount of gas (due to, e.g., ram or tidal 
stripping), it is clear that both assumptions are likely to 
hold until late times. Such galaxies may move along the 
correlation as they evolve, but they will remain on the 
same line. However, such galaxies are not the only ob- 
jects which may retain this correlation imtil the current 
cosmic epoch: in galaxies where star formation is halted 
early on (when our derivation assumptions are still valid) 
and gas is stripped from them, neither the metallicity nor 
the dark matter or the stellar mass of the galaxy would 
be affected. Hence the correlation between M/M* and 



Z would be preserved. This would then be a typical ex- 
ample of a "fossil" correlation, established in the early 
phases of the evolution of a galaxy and preserved until 
the present cosmic epoch due to lack of additional evo- 
lution of the observable quantities involved. 

Figure 6 shows the mass- weighted metallicity of stars of 
the simulated galaxies as a function of stellar mass. The 
correlation of the two quantities is very tight and can be 
described by a power law at small masses. At large stellar 
masses, M* > 5 x 10^ M©, the correlation becomes shal- 
lower. This is qualitatively consistent with the flattening 
of the Z — Msr relation found in the SDSS data (Tremonti 
et al. 2004) , although in their observations the break oc- 
curs at much larger masses: m ^ 3x 10^° Mq. Note, how- 
ever, that this figure should not be directly compared to 
the present-day correlation over a wide range of stellar 
masses because massive galaxies will undergo substan- 
tial evolution in both their stellar mass and metallicity 
between z = 3 and the present, while the evolution of 
dwarf galaxies is expected to be much slower. For some 
of the dwarf galaxies the evolution can be halted at high 
redshifts if their gas is removed by some process (e.g., 
ram pressiue stripping). Thus, the comparison with the 
correlations exhibited by the local galax:ies is most mean- 
ingful for smaller mass objects. 

The solid lines in Figure 6 show the results of our ana- 
lytic model. The model reproduces the slope of the corre- 
lation in the two mass regimes, the location of the break, 
and the observed redshift evolution seen in the simulated 
galaxies. The correlation in the model arises mainly due 
to the inefficiency of star formation in dwarf galaxies and 
the existence of star formation threshold, without gas 
outflows. Interestingly, direct extrapolation of the model 
to z = indicates that the flattening of the M, — Z re- 
lation shifts to M* ~ few x 10^° M©, in agreement with 
the SDSS data. However, the reader should keep in mind 
that certain assumptions {Z ^ Zq, ^ Mgas) of our 
analytic model are expected to eventually break down at 
late times (especially for the higher-mass objects), and 
hence additional effects may become important in shap- 
ing the present-day properties of higher-mass galaxies. 

5.3. Correlations Involving Stellar and Total Masses 

Figure 7 shows the stellar surface density, /x*, against 
galaxy stellar mass. The surface density is calculated 
from the simulation data as /it* oc M*/_R^, where i?* is 
the half-light radius, deflned as the radius which includes 
50% of the stellar mass. The data occupy a primary locus 
in the /i* — M* plane, which is traced by the prediction 
of our analytic model, in which stellar surface density 
is calculated as /z oc M^/r^j^. At masses M* < 10^ 
there are several outliers, as in this regime star formation 
becomes more stochastic. As for the other considered 
correlations, we see that the results of different simula- 
tions are consistent, although the NFEC and FNEC-RT 
simulations do show a higher /i* for comparable at 
the high-M* end, an indication that in these simulations 
the stellar populations appear to be more concentrated. 
We should note however that stellar surface densities are 
probably less reliable than global properties, such as stel- 
lar mass, because they depend on the internal structure 
of the stellar distribution in galax:ies and may be more 
susceptible to resolution effects. 

Figure 8 shows the stellar mass of the objects in our 
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Fig. 5. — Ratio between total mass and stellar mass (assumed to be proportional to the mass-to-light ratio of each object), plotted 
against stellar metallicity. Left panel: simulation FEC; middle panel: simulation NFEC; right panel: FNEC-RT simulation. Different 
symbols correspond to different redshifts as detailed in the legend. The solid line corresponds to the Mtot/Mt oc scaling 
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Fig. 6. — Stellar metallicity as a function of stellar mass. Left panel: simulation FEC; middle panel: simulation NFEC; right panel: 
FNEC-RT simulation. Different symbols correspond to different redshifts as detailed in the legend. The solid and dashed lines show the 
results of our analytic model for z = 4 and z = 9, respectively. 
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Fig. 7. — Stellar mass within ii« over (where R* is the half-light radius, which includes 50% of the stellar mass of each object), 
plotted against stellar mass. Left panel: simulation FEC; middle panel: simulation NFEC; right panel: FNEC-RT simulation. Different 
symbols correspond to different redshifts as detailed in the legend. The solid line corresponds to the prediction of out analytic model for 
z = 4. 



10 



simulations as a function of their maximum rotational ve- 
locity (the TuUy-Fisher relation). A break of the scaling 
law exists at a stellar mass of ~ 10^ Mq. For objects with 
stellar masses above the break, the stellar mass scales as 
with the best fit slopes of 3.41 ± 0.33, 3.34 ± 0.33, 
and 3.37 ±0.22 for the PEG, NFEC, and FNEC-RT sim- 
ulations. Below the break the Vm — Af* relation steepens 
considerably. The location of objects on the M, — Vm 
plane exhibits a considerable shift with redshift. This is 
a trivial consequence of the redshift-dependence of the 
virial density. 
Prom the virial theorem we have 

V^(x^/M/R~. (7) 

At the high redshifts under consideration, the universe 
is matter-dominated, and the virial overdensity d^iv = 
Pvir/Pm — 1 ~ IStt^ is csscntially independent of z. The 
virial radius is then determined from the relation 

^7ri?3.^p^(^)(l + 5vir) =M. (8) 

Given that Pmiz) = Pm,o(l + z)^, we get 

Ml/3 

Rvir CX —— . (9) 
1 + 2: 

and hence 

V^(xM^/^VT^. (10) 

In the lower panels of Figure 8, we have scale d the V m 
values of all data points of the upper panel by ^1 + z to 
correct for this redshift evolution. Using M oc Afgas and 
the Af* = Af*(A{fgas) dependence calculated using our 
analytic model, we get the solid and dashed lines shown 
in Figure 8, for z = 4 and z = 9, respectively. Once 
again, the model describes the results of the simulation 
very well. 

Finally. Figure 9 shows the baryonic TuUy-Fisher re- 
lation, where the total baryonic mass (stars and gas) 
of each object is plotted against Vm- In this case, the 
-^baryon OC scaliug (represented by the solid line) 
shows no evidence of a break. This is similar to the 
behavior of the observed dwarfs, which exhibit a break 
in the M^, — Vm relation, but a tight power law A^b — Kn 
relation (see § 2). 

5.4. The Effective Yield 
The effective yield, defined as 

^ WKY ^''^ 

has been widely used as a diagnostic of the evolution of 
the baryonic component of galaxies, and more specifi- 
cally as a "litmus test" of the validity of the closed-box 
approximation. Here, /g ~ A/gas/(Af* + A/fgas) is the 
fraction of baryons in the gas phase. Under the closed- 
box assumption, the effective yield is always equal to the 
true yield ytruei defined as the mass in newly synthesized 
metals returned to the ISM by a stellar population nor- 
malized to the stellar mass of this population locked-up 
in stellar remnants and long-lived stars. ^ 

^ The true yield is related to the mass of newly synthesized 
metals normalized to the total mass of the parent generation of 
stars, qz- J/true = Qz/o-, where a is the lock-up fraction. 



Edmunds (1990) showed that yea cannot exceed ytruo- 
Values of y^s lower than ytruo arc indicative of the devia- 
tions from the closed-box model due to either inflows or 
outflows of gas and metals. Observationally, dwarf galax- 
ies {Vm ^ lOOkm/s) have low values of j/eff- This has 
been interpreted as strong evidence for significant out- 
flows of metals and/or gas, which can reduce the metal- 
licity and/or the gas fraction and are expected to be more 
prevalent in smaller objects with shallower gravitational 
potentials (Garnett 2002; Tremonti et al. 2004; Pilyugin 
et al. 2004). 

In the previous section we showed that global correla- 
tions of galaxy properties arise in our simulations without 
SN-driven outflows. We have shown that these correla- 
tions can be reproduced with an open-box evolution with 
cosmological mass inflow but no outflows and 1) mass- 
dependent gas distribution given by equation (B3) and 2) 
star formation obeying the Kcnnicutt law with a surface 
density threshold. In this model the correlations, similar 
to those of observed galaxies, arise due to the increasing 
fraction of gas at surface densities below the star forma- 
tion threshold density (and, correspondingly, inefficiency 
of star formation) in smaller mass objects. This trend, 
in turn, is due to inefficiency of cooling and more ex- 
tended gas distributions for smaller masses. Given that 
the observed values of the effective yield are considered 
to be evidence for outflows, in this section we compare 
the predictions of our simulations for the values of the 
effective yield. 

Upper panels of Figure 10 show the effective yield as 
a function of circular velocity for the objects in simula- 
tions PEG (with feedback) and NPEG (no feedback) at 
z = 4. Here y^g? is calculated taking into account all 
of the gas and metals inside the virial radius of each ob- 
ject. As expected, in the absence of outflows the effective 
yield is constant and is close to the true yield adopted in 
simulations. 

However, we should recognize that in observations the 
gas fraction only includes gas of densities and tempera- 
tures accessible through HI observations (e.g., Tremonti 
et al. 2004; Garnett 2002; Pilyugin et al. 2004; Lee et al. 
2006a). If a fraction of metals resides in the gas phase 
that is not observed, the measured effective yield will 
be lower, just as in the case when the enriched gas is 
ejected from the system completely. The lower panels 
in Figure 10 show the values of calculated using 
metals only within the stellar extent, defined as the ra- 
dius that includes 90% of the total stellar mass, and cold 
gas only (T < lO^'K). The metallicity entering the "ob- 
servable" effective yield is calculated as the ratio of the 
metal mass within the 90% stellar mass radius, and the 
gas mass within the same radius. The gas entering the 
gas fraction is all gas within the halo with T < lO^K. In 
this case, there is a marked decrease of with decreas- 
ing circular velocity at < 100 km s"^. The values of 
the effective yield for the small-mass systems are now 
comparable to those measured for dwarf galaxies. The 
decreasing effective yield for lower mass systems suggests 
that a fraction of metals in these dwarf galaxies is not 
in the observed cold gas. Given that the effective yield 
within the virial radius is constant and is equal to its 
true value, the extra metals must reside within the virial 
radius but in a hot phase gas, detectable perhaps only 
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Fig. 8. — Maximum circular velocity as a function of stellar mass (Tully-Fisher relation). Left panel: simulation FEC; middle panel: 
simulation NFEC; right panel: FNEC-RT simulation. Different symbols correspond to different redshifts as detailed in the legend. The 
solid and dashed lines show the results of our analytic model for z = 4 and z = 9, respectively. In the lower row plots we have shifted Vm 
by a factor of (1 + z)~^-^ to correct for the explicit redshift dependence entering Vm through its dependence on the virial radius. 



by its absorption in spectra of background sources and 
which is thus not included in the calculation of the gas 

fraction. 

In order to reproduce the correct trend with galaxy 
mass, low-mass objects should have a larger fraction of 
the metals residing in the unobservable gas phase. Figure 
1 1 shows the fraction of the metal mass within the cold 
gas disk, f^°^'^ and shows that our simulated galaxies do 
indeed exhibit such a trend. For higher-mass galaxies, 
this fraction is close to unity, indicating that most of the 
metals are in the cold gas phase. For low-mass galax- 
ies, however, this fraction is significantly below unity, 
which means that a significant fraction of metals is in 
the hotter, unobservable gas. It is this trend that is 
responsible for the decreasing effective yield values for 
smaller mass galaxies in the bottom panels of Figure 10. 
The left and right columns, corresponding to simulations 
FEC and NFEC respectively, show similar trends, which 
shows that supernova-driven outflows are not responsi- 
ble in our simulations, and therefore are not necessarily 
responsible in nature, for the presence of metals outside 
of the cold phase. 

There is widespread turbulence in the interstellar 
medium of our simulated galaxies, driven by gravita- 
tional instabilities, cold accretion flows (Birnboim & 
Dekel 2003; Keres et al. 2005) reaching well inside the 
virial radius and stirring up the gas. However, the domi- 
nant source of both turbulence in the cold gas as well as 
mixing of metals in the hot phase are most likely merg- 



ers with other galaxies and dark matter subhalos, fre- 
quent at high redshifts. There is also the possibility that 
numerical diffusion also contributes to some extent to 
the efficient spread of metals within the virial radius of 
each object. Higher resolution simulations are required 
to quantify this effect. 

If this picture is also applicable in nature, the pre- 
diction of our model would be that, in the case of small 
dwarfs with low values of the observed effective yield, the 
gas well outside of the extent of the observable cold disk 
should be significantly enriched with heavy elements. Al- 
though to our knowledge there are no current observa- 
tional constraints on the metallicity of such gas^°, metal- 
licity measurements in absorption and tests of our con- 
jecture may be possible in the future. The current mea- 
surements for ~ galaxies indicate that cold gas is sig- 
nificantly enriched to large radii where little or none in 
situ star formation should be occuring (Chen et al. 2005). 
In addition, direct measurements of metal abundances at 
different radii show that metallicity gradients in dwarfs 
galaxies are weak to at least four exponential disk scale 

Some existing measurements of the metallicity of neutral HI 
gas for dwarf galaxy I Zw 18 (e.g., Aloisi et al. 2003; Lccavclier 
des Etangs et al. 2004) indicate that the metallicity of the neutral 
gas is comparable or somewhat below that of the HII gas for which 
the metallicities are typically measured. These measurements of 
the neutral phase metallicity, however, are still limited to regions 
within the stellar extent. It is not yet clear whether results for 
this very low metallicity galaxy are also true for dwarf galaxies in 
general. 
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Fig. 9. — Maximum circular velocity as a function of baryonic mass {baryonic TuUy-Fisher relation). Left panel: simulation FEC; middle 
panel: simulation NFEC; right panel: FNEC-RT simulation. Different symbols correspond to different redshifts as detailed in the legend. 
The solid line corresponds to the Mj,a^yojj oc scaling. In the lower row plots we have shifted Vm by a factor of (1 + z)""'® to correct for 
the explicit redshift dependence entering Vm through its dependence on the virial radius. 



lengths (Lee et al. 2006b, and references therein). This 
may indicate mixing of metals in the gaseous disks is 
indeed efficient. 

6. DISCUSSION 

Before a detailed comparison with observational results 
is attempted, it is important to consider how the results 
that we obtain for simulated galaxies at high-redshifts 
can relate to the galaxies observed at lower redshifts. Our 
analytic model shows that the correlations we observe in 
the simulated galaxies arise because galaxies are gas rich 
and star formation is increasingly inefficient for lower 
mass systems (which is due to the existence of density 
threshold for star formation). These results thus are not 
applicable to galaxies which formed significant fraction 
of their stars under different conditions. 

However, we argue that our results are applicable to 
certain types of the nearby dwarf galaxies. All of such 
galaxies could have undergone outflow-free evolution in 
the gas-rich regime at high redshift, similar to the our 
simulated galaxies. Their subsequent fates could be 
roughly classified into three scenario. Some galaxies 
could have continued to grow their mass and form stars 
and c:onsumed most of their gas. Such galaxies could 
have formed most of their stars in the regime where our 
results are not applicable. This scenario likely applies for 
those small-mass, high-2 systems that have experienced 
significant mass growth and became massive objects to- 
day. Another scenario is that star formation could have 



been truncated at high redshifts, while still in the gas- 
rich regime, by ram pressure and tidal forces from the 
nearby massive galaxy. Such scenario is probably appli- 
cable for the nearby dwarf spheroidal galaxies (Kravtsov 
et al. 2004; Mayer et al. 2007). Finally, the third scenario 
is that the mass of a high-z dwarf system remained low 
throughout its evolution. In this case most of the stars 
could have formed in the low-efficient mode of star for- 
mation and the galaxy may have remained in the gas-rich 
regime until z = 0. 

In either one of the latter two scenarios, high-redshift 
correlations between metallicity, total mass, stellar mass 
and surface brightness would be preserved to 2; = 0. Note 
that both nearby gas-poor dwarf spheroidals and gas- 
rich dwarf irregulars exhibit similar correlations of their 
global properties such as luminosity and metallicity (e.g., 
Grebel et al. 2003). The slopes of these correlations are 
indeed very close to the ones we find in our models of 
high-redshift objects. This may indicate that the evolu- 
tion of these galaxies has proceeded according to the last 
two scenarios. 

It is thus interesting to examine how our simulation re- 
sults high-redshift results compare with observations of 
dwarf galaxies in the present-day universe. In Figure 12 
we plot observed properties of dwarf galaxies, the z = 4 
results from the FEC simulation, and the predictions of 
the analytic model for the same redshift. The upper 
left panel shows the metallicity as a function of stellar 
mass. Observational data are taken from Dekel & Woo 
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(2003) (in turn taken from Mateo 1998, van den Bosch 
2000 and references therein) and from (Lec ct al. 2006a, 
their Table 1). In the case of the Lee et al. (2006a) data, 
we used the conversion relation fitted by (Simon et al. 
2006, their Eq. 16) to convert 12 + log(0/H) to [Fe/H]. 
The upper right panel shows the surface brightness as a 
function of stellar mass. In the lower left panel, we plot 
the stellar mass as a function of the maximum circu- 
lar velocity (the TuUy-Fisher relation) , corrected against 
explicit evolution with redshift. Finally, in the lower 
right panel we show the effective yield as a function 
of maximum circular velocity. The overplottcd obser- 
vational points are taken from (Lee et al. 2006a, their 
Table 5), (Garnett 2002, their Table 4) and (Pilyugin 
et al. 2004, we have used data on oxygen abundances and 
gas fractions from their Table 7 to compute the effective 
yield). The effective yield in the simulated objects has 
been calculated using the metals and gas within the stel- 
lar disk to calculate the mctallicity, but all of the cold 
(< 10** K), Hl-observable gas in each halo to calculate 
/gas, so as to be directly comparable with the observa- 
tional data. Although the scatter is appreciable both 
within the same dataset but most prominently between 
different datasets^^, the correlations present in the sim- 
ulation results are generally very similar with the ones 
found in the observational data. 

We have shown in §5 that in all of our models, tight cor- 
relations, similar to those observed in present-day dwarfs, 
are developed already at very high redshifts between ob- 
servables of low-mass galaxies, and they are preserved 
at least until 2; ~ 3. Energy feedback for supernovae 
and galactic winds was not found to be either necessary 
or important for establishing these correlations. Our in- 
terpretation of these correlations in our models and the 
corresponding ones in observed systems is that they oc- 
cur as a result of the inefficiency of star formation in 
low-mass systems, due to inefficient cooling and the as- 
sociated low gas surface-densities (a concept consistent 
with observations suggesting that in low-mass galaxies 
the surface density of cold gas is typically low, e.g., van 
Zee et al. 1997; Martin & Kennicutt 2001; Auld et al. 
2006). As a result, stellar mass and metal enrichment 
are increasingly suppressed at lower masses, resulting in 
the observed correlations between stellar mass and other 
observables. The inefficiency of star formation below the 
critical gas surface density also helps to explain the ob- 
served gas fractions as a function of galaxy luminosity 
and surface brightness (van den Bosch 2000), the stel- 
lar truncation radii and gas extent in the observed disks 
(van den Bosch 2001), the shallow slope of the faint-end 
of the galaxy luminosity function (Verde et al. 2002; Cro- 
ton et al. 2006), the paucity of faint dwarf galaxies in the 
Local Group (Verde et al. 2002; Kravtsov et al. 2004), 
and observations of extreme Hl-dominated dwarfs (Be- 
gum et al. 2005, 2006; Warren et al. 2004, 2006), which 
are not deficient in baryons but are nevertheless underlu- 
minous. These results imply that the star formation law 
governing conversion of cold gas into stars, rather than 
SN-driven outflows, is the dominant factor in shaping 

Pilyugin ct al. (2004) and Garnett (2002) have a large fraction 
of their samples overlapping, often deriving very different rnetal- 
licities for individual cases. Thus, the scatter between different 
datasets in the case of the effective yield is representative of sys- 
tematic uncertainties in such determinations. 



properties of faint galaxies. 

In addition, there are observational indications that 
feedback does not in fact drive significant amounts of gas 
out of low-mass galaxies. Given that lower-mass galax- 
ies would be more susceptible to strong supernova-driven 
gas outfiows, one could expect that gas fractions should 
decrease with decreasing stellar mass. However, the op- 
posite trend is observed (McGaugh & de Blok 1997; Bell 
& de Jong 2000; Garnett 2002; Geha et al. 2006). In 
fact, many of the smallest galaxies have gas fractions in 
excess of 0.8-0.9 (Geha et al. 2006; Lee et al. 2006a). 

At the same time however, there is indirect evidence for 
the importance of metal loss in winds. The high mctal- 
licity of diffuse intcrgalactic gas in groups and clusters 
(e.g., Renzini et al. 1993) and in the intergalactic medium 
around galaxies (e.g., Adelberger et al. 2005) is gener- 
ally attributed to the action of winds. The winds may 
also play a role in faster evolution of the stellar mass- 
metallicity relation (De Rossi et al. 2006). 

An interesting additional constraint comes from the 
observed HI mass function (HIMF) of galaxies (Zwaan 
et al. 2005). Its small-mass end is significantly shallower 
than the expected mass function of dark matter halos, 
implying that some sort of suppression of the gas con- 
tent in halos is needed. At this point, there is no ac- 
cepted explanation of the shallow small-mass end of the 
HIMF, although Mo et al. (2005) have recently argued 
that it can be explained by a combination of gas heating 
by the cosmic UV background and preheating of accret- 
ing gas by filaments and pancakes in which galaxies are 
embedded. It is interesting that regardless of the nature 
of suppression mechanism it should be almost indepen- 
dent of halo mass. Strong mass dependence of gas sup- 
pression would likely introduce a break in the baryonic 
TuUy-Fisher relation, which is not observed. 

Prom a chemical evolution point of view, outflows 
would constitute a departure from closed-box evolution. 
The qtiantity that has been traditionally used to quan- 
tify the mctallicity evolution of a galaxy and its devia- 
tion from the closed-box model is the effective yield, yefi 
(Garnett 2002; Pilyugin et al. 2004; Tremonti et al. 2004). 
Lower-mass galaxies (Vm ^ lOOkm/s) have in fact been 
observed to have an effective yield lower than the closed- 
box model prediction for their observed gas fraction, and 
this has been traditionally interpreted as an indication 
for metal outflows. Recently, Dalcanton (2006) showed 
that outflows of gas with the same average metallicity as 
the host galaxy can only have a minor effect on the effec- 
tive yield and cannot reproduce the observed reduction of 
UeB with respect to the closed-box value in low-mass ob- 
jects, and suggested that metal-enhanced outflows may 
instead be responsible for this effect. 

In this paper, we show that there is an alternative 
explanation for the trend for decreased effective yields 
at the lowest-mass objects. The effective mixing in the 
ISM at these high redshifts spreads the metals produced 
through star formation to larger regions in small galaxies 
than are accessible through observations. When the ef- 
fective yield in our simulated galaxies is calculated taking 
into account gas, stars and metals out to the virial ra- 
dius of each object, no dependence of the effective yield 
with mass is observed. Conversely, when we calculate 
the effective yield taking into account only the gas, stars 
and metals that would be accessible to observations, a 
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trend of decreasing effective yield with mass is revealed, 
similar to the one seen in observational data. This trend 
becomes much less pronounced in higher-mass galaxies, 
consistent with the findings of Erb ct al. (2006) who 
found no appreciable dependence of effective yield on 
stellar mass for objects with M* > lO^M© and for 2; ~ 2. 

We would like to stress again that gas outflows escap- 
ing the host halo are not required to explain this trend, 
although they may affect the exact form of the relation 
(normalization and slope) and its evolution. Note how- 
ever that gas flows within the host halo may be important 
in spreading the produced metals to regions considerably 
more extended than the stellar disk. The consistency be- 
tween results from our simulations with and without en- 
ergy feedback from supernovae suggests that some other 
process drives turbulence and large-scale motions of the 
interstellar medium, such as gravitational instabilities, 
galaxy mergers, or cold accretion flows. This explana- 
tion of the trend of effective yield with galaxy mass can 
be potentially tested by measurements of the metallicity 
of the gas in the outskirts of the gaseous disks of dwarf 
galaxies. 

7. CONCLUSIONS 

In this paper, we use a suite of high-resolution cos- 
mological simulations to examine the role of different 
physical processes in establishing the observed scaling 



relations of dwarf galaxies. All simulations included a 
recipe for star formation and metal enrichment of the 
ISM. Three simulations (FEC, FNEC-RT, F2NEC-RT) 
additionally included energy feedback from supernovae 
to the ISM. Simulations FNEC-RT and F2NEC-RT fol- 
lowed detailed 3-D radiative transfer from individual UV 
sources and the chemistry network of ionic species of hy- 
drogen, helium, and molecular hydrogen to determine 
the radiative heating and cooling of the baryonic gas, 
while the other two simulations used tabulated heating 
and cooling rates and a uniform ionizing background. In 
order to interpret our simulation results, we develop and 
use an analytic model, which follows the evolution of the 
dynamical mass, gas, and stellar components as well as 
the metallicity of galaxies. We assumed an open box 
model with no outflows but with accretion of mass at 
a rate proportional to the galaxy mass, and star forma- 
tion following the Kennicutt law with a critical density 
threshold for star formation. Our main results and con- 
clusions can be summarized as follows. 

• In all simulations, correlations between global 
quantities, such as the stellar mass-metallicity rela- 
tion, very similar to those observed for the nearby 
dwarf galaxies arise in simulated galaxies with stel- 
lar masses lO^M© < < lO^M© for z > 3.3. 

• In galaxies of higher masses these correlations ex- 
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hibit a break and the dependence of global proper- 
ties on stellar mass flattens off. A similar flatten- 
ing of the correlations is also found in observational 
data, although the location of the break may have 
evolved to higher masses by the present epoch (e.g., 
Tremonti et al. 2004). 

• The key finding of this study is that neither the 
inclusion of supernova energy feedback, nor the in- 
clusion of 3-D radiative transfer significantly affects 
any of these correlations, a strong indication that 
thermal and radiative processes (as well as associ- 
ated galactic outflows and winds) are not required 
for the observed correlations to arise. Using our 
analytic model, we show that these correlations in 
the simulated galaxies arise due the increasing in- 
efficiency of star formation with decreasing mass 
of the object. The inefficiency of conversion of gas 
into stars and its dependence on mass are due to 
the critical gas column density threshold of the star 
formation law. 

• We thus argue that trends similar to those exhib- 
ited by observed dwarf galaxies can arise without 
loss of gas and metals in winds, although winds 
may affect the exact form of the relation (normal- 
ization and slope) and its evolution. 

• The observed trends of a decreasing effective yield 
with decreasing galaxy mass can be reproduced in 

our simulated galaxies reasonably well, and with- 
out winds, if the effective yield is calculated taking 
into account only the gas, stars, and metals that 
would be accessible to observations. We show that 
this trend is due to efficient mixing in the ISM, 
driven by disk instabilities, mergers, and possibly 
cold accretion flows. The efficient mixing redis- 



tributes the metals far beyond the stellar extent in 
smaller mass galaxies. Thus a significant fraction 
of metals in small galaxies is outside the stellar ex- 
tent leading to a decrease of the estimated yield 
if only metals within the stellar radius are taken 
into account. This explanation can potentially be 
tested by measurements of the gas metallicity in 
the outskirts of the gaseous disks of dwarf galaxies. 
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APPENDIX 
GAS MASS ACCRETION RATE 

In the spherical collapse model, an overdensity of mass M at a time t accretes matter at a rate F = dM/dt cx M 
(Fillmore & Goldreich 1984; Bertschinger 1985).^^ Such a relation is also consistent with the results of cosmological 
simulations, in which halo mass grows with time as M{z) oc exp[— C2:] (Wechsler et al. 2002), which gives dM/dt = 
f(z)M. In general, f{z) depends nontrivially on cosmic time (with the exception of the Einstein-de Sitter miiverse 
where the evolution is self-similar and the / scales with time simply as cx i^^), as well as the cosmological parameters. 

Since we consider early epochs in the evolution of galaxies we assume that the baryon mass is dominated by gas. 
For example, for the simulated halos in the redshifts of interest: Mgas > lOM*. We assume thus that Mgas ~ Mb = 

ObM/l^m, so 

F = f{z)M^ . (Al) 
An expression for f{z) is obtained by using the following simple recipe: 

F ^ 47r[rvir(z)]^(l + <5cxt)/0m,z(z)Wacc(z) , (A2) 

where Pm,z — Pm,o(l + z)'^ if^ the mean matter density of the universe at redshift z, rvir(z) = [3A'//47r/9ni.2(l + (^vir)]"'^^^ 
is the virial radius of an object of mass M at redshift z, Jvir is the virial overdensity (equal to IStt^ w 178 for the 
high redshifts considered here), 1 4- (5cxt is the local compression factor right outside the virial radius of the object, 
and Wacc(-z) — ■\/2GM/rvir(z) is the free-fall velocity Equation (A2) is only approximate, and we use the value of 
(1 + (^cxt) to set its exact normalization so that it adequately reproduces our simulations results, which we find to be 
the case for 1 + 5ext ~ 20. Similar values of the external compression factor used in analytic cosmological accretion 
models yield results in general agreement with the accretion rates and energetics of objects in cosmological simulations 
(Pavhdou & Fields 2006). 

Although analytic solutions were derived for an Einstein - de Sitter universe, the result holds under spherical collapse in any cosmological 
model. 

13 Note that r-vir <x M^/^ and Uacc oc M'^/'^r'-^J'^ oc mV3 so F oc r^j^uacc oc M. 
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STAR FORMATION RATE 

To calculate '0(Mgas) we assume that gas in each halo settles into an exponential disk obeying the empirical Kennicutt 

star formation law^^, with stars forming only above a certain surface gas density threshold. Specifically, we assume 
that the star formation rate per surface area, S*, depends on the surface density of gas, Sg, as 



2.5 X 10" 



1.4 



.Mq pC-2 

(Kennicutt 1998). The surface gas density profile is assumed to be 

Sg(r) = Soexp 



Mgkpc Vr"^ 



(Bl) 



(B2) 



Following Kravtsov et al. (2004), we adopt the following scaling of the disk scale radius, ra, with the total halo mass: 

2/3 

(B3) 



rd{z) = 2 ^/^Arvir(-z)exp 



where M4 = 2.2 x IQ^^Mq is the total mass of a halo with virial temperature IQ^K virializing at 2; = 0, and A is the 
angular momentum parameter, which has a log-normal distribution of values with scale and shape parameters 0.045 
and 0.56 respectively (e.g., Vitvitska et al. 2002, with the peak of the distribution at A = 0.033). The finite spread in 

possible angular momentum parameters contributes to the scatter in the observed correlations. We take c = 0.08 by 
calibrating against our simulations data. The rcdshift dependence in the exponential enters through the dependence 
of the virial temperature on the redshift of virialization. 

The scaling given by equation (B3) adopts the commonly used scaling (r^ oc Aryir) for massive halos (Mo et al. 
1998), while for systems with the virial temperatures Tvir ^ lO^K, it assumes that the gas cannot cool efficiently and 
cannot settle into a rotationally supported disk but settles into a more extended equilibrium configuration in the halo 
(see Kravtsov et al. 2004). As a result, the associated gas surface density decreases steeply with decreasing mass and 
star formation is inefficient or is completely suppressed in smaller objects. It is this ineflaciency of star formation in 
dwarf galaxies that is primarily responsible for the observed correlations. 

Star formation is suppressed for gas surface densities Eg > ~ SMqpc^^ (Martin & Kennicutt 2001), which 
implies the threshold radius of 

rth(^)=rd(2)ln^. (B4) 

This radius decreases with decreasing galaxy mass. In other words, star formation proceeds within smaller radii in 
galaxies of smaller masses. 
The star formation rate V' can then be calculated as follows: 



V'(Mgas,2)= / S*27rrdr 
Jo 



= 8 X lO-^'Mgyr" 



1 ( -Tdiz) 
kpc 



1 - I 1 + — ^ exp 



So(^) 

Mq pC-2 

_L4rth(£) 
rd{z) 



1.4 



(B5) 



In turn, Sq and its redshift dependence can be obtained through 

poo 

Mgas= / 27rrrfrSg(r) = 27r[rd(2)]^So(2) . (B6) 

Jo 

where is obtained as a function of Mgas from equation (B3) assuming Mgas ~ flb-^/f^m- 

Using equation (Al) for the accretion rate and equation (B5) for the star formation rate, we can integrate equations 
(3) and (4) numerically. As we show in §5, the agreement of the resulting dependence of M* on Mgas with the relation 
obeyed by the simulated galaxies (see Fig. 3) is striking, especially when one takes into account that the values of the 
parameters A and Sth are either on or very close to their "most probable" values. 
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